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, Abstract 

We study the PageRank mass of principal components in a bow-tie Web Graph, as 
(f) ■ a function of the damping factor c. Using a singular perturbation approach, we show 

that the PageRank share of IN and SCC components remains high even for very large 
values of the damping factor, in spite of the fact that it drops to zero when c — » 1. 
However, a detailed study of the OUT component reveals the presence "dead-ends" 
(small groups of pages linking only to each other) that receive an unfairly high ranking 
when c is close to one. We argue that this problem can be mitigated by choosing c as 
small as 1/2. 

1 Introduction 

> , 

The link-based ranking schemes such as PageRank pQ, HITS [2J, and SALSA [3] have been 
successfully used in search engines to provide adequate importance measures for Web pages. 
^ In the present work we restrict ourselves to the analysis of the PageRank criterion and use 

the following definition of PageRank from 0]. Denote by n the total number of pages on 
the Web and define the n x n hyper-link matrix W as follows: 



o 
o 



- 1—^ 

X 




if page i links to j, 
if page i is dangling, 
otherwise, 



for i,j — l,...,n, where di is the number of outgoing links from page i. A page is called 
dangling if it does not have outgoing links. The PageRank is defined as a stationary distri- 
bution of a Markov chain whose state space is the set of all Web pages, and the transition 
matrix is 

G = cW + (1 -c)(l/n)l T l. (2) 

Here and throughout the paper we use the symbol 1 for a column vector of ones having by 
default an appropriate dimension. In ([2]), 1 T 1 is a matrix whose all entries are equal to one, 
and c £ (0, 1) is the parameter known as a damping factor. Let tt be the PageRank vector. 
Then by definition, irG — tt, and ||7r|| = ttI = 1, where we write ||x|| for the Li-norm of 
vector x. 
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The damping factor c is a crucial parameter in the PageRank definition. It regulates 
the level of the uniform noise introduced to the system. Based on the publicly available 
information Google originally used c = 0.85, which appears to be a reasonable compromise 
between the true reflection of the Web structure and numerical efficiency (see [5] for more 
detail). However, it was mentioned in [6] that the value of c too close to one results into 
distorted ranking of important pages. This phenomenon was also independently observed 
in [7j. Moreover, with smaller c, the PageRank is more rebust, that is, one can bound the 
influence of outgoing links of a page (or a small group of pages) on the PageRank of other 
groups [3] and on its own PageRank [TJ. 

In this paper we explore the idea of relating the choice of c to specific properties of 
the Web structure. In papers [9] [10] the authors have shown that the Web graph can be 
divided into three principle components. The Giant Strongly Connected Component (SCC) 
contains a large group of pages all having a hyper-link path to each other. The pages in 
the IN (OUT) component have a path to (from) the SCC, but not back. Furthermore, the 
SCC component is larger than the second largest strongly connected component by several 
orders of magnitude. 

In Section [3] we consider a Markov walk governed by the hyperlink matrix W and ex- 
plicitly describe the limiting behavior of the PageRank vector as c — > 1. We experimentally 
study the OUT component in more detail to discover a so-called Pure OUT component 
(the OUT component without dangling nodes and their predecessors) and show that Pure 
OUT contains a number of small sub-SCC's, or dead-ends, that absorb the total PageRank 
mass when c = 1. In Section |4] we apply the singular perturbation theory [IT] , [T2 ] [T3 ] [14] 
to analyze the shape of the PageRank of IN+SCC as a function of c. The dangling nodes 
turn out to play an unexpectedly important role in the qualitative behavior of this func- 
tion. Our analytical and experimental results suggest that the PageRank mass of IN+SCC 
is sustained on a high level for quite large values of c, in spite of the fact that it drops to 
zero as c — > 1. Further, in Section [5] we show that the total PageRank mass of Pure OUT 
component increases with c. We argue that c = 0.85 results in an inadequately high ranking 
for Pure OUT pages and we present an argument for choosing c as small as 1/2. We confirm 
our theoretical argument by experiments with log files. We would like to mention that the 
value c = 1/2 was also used in [15] to find gems in scientific citations. This choice was 
justified intuitively by stating that researchers may check references in cited papers but on 
average they hardly go deeper than two levels. Nowadays, when search engines work really 
fast, this argument also applies to Web search. Indeed, it is easier for the user to refine a 
query and receive a proper page in fraction of seconds than to look for this page by clicking 
on hyper-links. Therefore, we may assume that a surfer searching for a page, on average, 
does not go deeper than two clicks. 

The body of the paper contains main ideas and results. The necessary information from 
the perturbation theory and the proofs are given in Appendix. 

2 Datasets 

We have collected two Web graphs, which we denote by INRIA and FrMathlnfo. The Web 
graph INRIA was taken from the site of INRIA, the French Research Institute of Informatics 
and Automatics. The seed for the INRIA collection was Web page www.inria.fr. It is a 
typical large Web site with around 300.000 pages and 2 millions hyper-links. We have 
collected all pages belonging to INRIA. The Web graph FrMathlnfo was crawled with the 
initial seeds of 50 mathematics and informatics laboratories of France, taken from Google 
Directory. The crawl was executed by Breadth First Search of depth 6. The FrMathlnfo 
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Web graph contains around 700.000 pages and 8 millions hyper-links. Because of the fractal 
structure of the Web [16] we expect our datasets to be enough representative. 

The link structure of the two Web graphs is stored in Oracle database. We could store 
the adjacency lists in RAM to speed up the computation of PageRank and other quantities 
of interest. This enables us to make more iterations, which is extremely important when 
the damping factor c is close to one. Our PageRank computation program consumes about 
one hour to make 500 iterations for the FrMathlnfo dataset and about half an hour for 
the INRIA dataset for the same number of iterations. Our algorithms for discovering the 
structures of the Web graph are based on Breadth First Search and Depth First Search 
methods, which are linear in the sum of number of nodes and links. 



3 The structure of the hyper-link transition matrix 

With the bow-tie Web structure [9] [TU] in mind, we would like to analyze a stationary 
distribution of a Markov random walk governed by the hyper-link transition matrix W 
given by |T]). Such random walk follows an outgoing link chosen uniformly at random, and 
dangling nodes are assumed to have links to all pages in the Web. We note that the methods 
presented below can be easily extended to the case of personalized PageRank [17] , when after 
a visit to a dangling node, the next page is sampled from some prescribed distribution. 

Obviously, the graph induced by W has a much higher connectivity than the original 
Web graph. In particular, if the random walk can move from a dangling node to an arbitrary 
node with the uniform distribution, then the Giant SCC component increases further in size. 
We refer to this new strongly connected component as the Extended Strongly Connected 
Component (ESCC). Due to the artificial links from the dangling nodes, the SCC component 
and IN component are now inter-connected and are parts of the ESCC. Furthermore, if 
there are dangling nodes in the OUT component, then these nodes together with all their 
predecessors become a part of the ESCC. 

In the mini-example in Figure [lj node represents the IN component, nodes from 1 to 
3 form the SCC component, and the rest of the nodes, nodes from 4 to 11, are in the OUT 
component. Node 5 is a dangling node, thus, artificial links go from the dangling node 5 to 
all other nodes. After addition of the artificial links, all nodes from to 5 form the ESCC. 




# 


INRIA 


FrMathlnfo 


total nodes 


318585 


764119 


nodes in SCC 


154142 


333175 


nodes in IN 








nodes in OUT 


164443 


430944 


nodes in ESCC 


300682 


760016 


nodes in Pure OUT 


17903 


4103 


SCCs in OUT 


1148 


1382 


SCCs in Pure Out 


631 


379 



Figure 1 : Example of a graph 



Figure 2: Component sizes in INRIA 
and FrMathlnfo datasets 



In the Markov chain induced by the matrix W, all states from ESCC are transient, that 
is, with probability 1, the Markov chain eventually leaves this set of states and never returns 
back. The stationary probability of all these states is zero. The part of the OUT component 
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without dangling nodes and their predecessors forms a block that we refer to as a Pure OUT 
component. In Figure [T] the Pure OUT component consists of nodes from 6 to 11. Typically, 
the Pure OUT component is much smaller than the Extended SCO However, this is the set 
where the total stationary probability mass is concentrated. The sizes of all components for 
our two datasets are given in Figure [2] Here the size of the IN components is zero because 
in the Web crawl we used the Breadth First Search method and we started from important 
pages in the Giant SCO For the purposes of the present research it does not make any 
difference since we always consider IN and SCC together. 

Let us now analyze the structure of the Pure OUT component in more detail. It turns 
out that inside Pure OUT there are many disjoint strongly connected components. All states 
in these sub-SCC's (or, "dead-ends") are recurrent, that is, the Markov chain started from 
any of these states always returns back to it. In particular, we have observed that there are 
many dead-ends of size 2 and 3. The Pure OUT component also contains transient states 
that eventually bring the random walk into one of the dead-ends. For simplicity, we add 
these states to the giant transient ESCC component. 

Now, by appropriate renumbering of the states, we can refine the matrix W by subdi- 
viding all states into one giant transient block and a number of small recurrent blocks as 
follows: 



W 



Qi 

Q m 

Rl • • • Rm T 



dead-end (recurrent) 

(3) 

dead-end (recurrent) 

ESCC+ [transient states in Pure OUT] (transient) 



Here for i — 1, . . . , to, a block Qi corresponds to transitions inside the i-th recurrent block, 
and a block Ri contains transition probabilities from transient states to the i-th recurrent 
block. Block T corresponds to transitions between the transient states. For instance, in 
example of the graph from Figure [TJ the nodes 8 and 9 correspond to block Qi, nodes 10 
and 11 correspond to block Q2, and all other nodes belong to block T. 

We would like to emphasis that the recurrent blocks here are really small, constituting 
altogether about 5% for INRIA and about 0.5% for FrMathlnfo. We believe that for larger 
data sets, this percentage will be even less. By far most important part of the pages is 
contained in the ESCC, which constitutes the major part of the giant transient block. 

Next, we note that if c < 1, then all states in the Markov chain induced by the Google 
matrix G are recurrent, which automatically implies that they all have positive stationary 
probabilities. However, if c = 1, the majority of pages turn into transient states with sta- 
tionary probability zero. Hence, the random walk governed by the Google transition matrix 
([2]) is in fact a singularly perturbed Markov chain. Informally, by singular perturbation we 
mean relatively small changes in elements of the matrix, that lead to altered connectivity 
and stationary behavior of the chain. Using the results of the singular perturbation theory 
(see e.g., [TTJ [121 EE HI]); m the next proposition we characterize explicitly the limiting 
PageRank vector as c — > 1 (see Appendix IA. 21 for the proof). 

Proposition 1 Let ~k vt,% be a stationary distribution of the Markov chain governed by Qi, 
2 = 1,..., to . Then, we have 



where 



lim7r(c) = [ir OUTl ■■■ ir OUTym 0] 

c— >1 



# nodes in block Qi 1 T . rrn—i- 



n n 



+ -l I [I-T]- 1 R i l)n 0UTti (4) 
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Figure 3: The PageRank mass of IN+SCC as a function of c. 



for i — l,...,m, I is the identity matrix, and is a row vector of zeros that correspond to 
stationary probabilities of the states in the transient block. 

The second term inside the brackets in formula ([4]) corresponds to the PageRank mass 
received by a dead-end from the Extended SCC. If c is close to one, then this contribution 
can outweight by far the fair share of the PageRank, whereas the PageRank mass of the 
giant transient block decreases to zero. How large is the neighborhood of one where the 
ranking is skewed towards the Pure OUT? Is the value c = 0.85 already too large? We will 
address these questions in the remainder of the paper. In the next section we analyze the 
PageRank mass IN+SCC component, which is an important part of the transient block. 



4 PageRank mass of IN+SCC 

In Figure 3 we depict the PageRank mass of the giant component IN+SCC, as a function of 
the damping factor, for FrMathlnfo. Here we see a typical behavior also observed for several 
pages in the mini- web from [6]: the PageRank first grows with c and then decreases to zero. 
In our case, the PageRank mass of IN+SCC drops drastically starting from some value c 
close to one. We can explain this phenomenon by highlighting the role of the dangling nodes. 

We start the analysis by subdividing the Web graph sample into three subsets of nodes: 
IN+SCC, OUT, and the set of dangling nodes DN. We assume that no dangling node 
originates from OUT. This simplifies the derivation but does not change our conclusions. 
Then the Web hyper-link matrix W in ([!]) can be written in the form 



W 



Q 
R P S 
Ill T ±11 T ill T 



OUT 

IN+SCC 

DN 



where the block Q corresponds to the hyper-links inside the OUT component, the block 
R corresponds to the hyper-links from IN+SCC to OUT, the block P corresponds to the 
hyper-links inside the IN+SCC component, and the block S corresponds to the hyper-links 
from SCC to dangling nodes. In the above, n is the total number of pages in the Web graph 
sample, and the blocks 11 T are the matrices of ones adjusted to appropriate dimensions. 

Dividing the PageRank vector in segments corresponding to the blocks OUT, IN+SCC 
and DN, 



' OUT " IN + SCC "DNjj 
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we can rewrite the well-known formula (see e.g. |18| ) 

n= LLE 1 T [I _ cP] -i (5) 

as a system of three linear equations: 

ttout[^ - cQ] ~ 7i"iN + sc C ci? - -tt dn 11 t = - — -1 T , (6) 

n n 

K m+S cc[I - cP] - -tt dn 11 t = — 1 T , (7) 

n n 

- TTlN+SCcCS' + 7T DN 7T DN 11 = 1 . (8) 

n n 

Solving (JH] [HD f° r ttin+scc we obtain 



(1 — c)a |~ T „ c 2 a 



1 — cp 



I-CP-- -Sill* 

1 — cp 



(9) 



where 

a = l/iV + SCCI/n and (3 = \DN\/n 

are the fractions of nodes in IN+SCC and DN, respectively, and u IN+scc = |/A r + S'CC|~ 1 1 T 
is a uniform probability row- vector of dimension \IN + SCC\. The detailed derivation of 
([9]) can be found in Appendix IA. 21 
Now, define 

(1 — Cid ccx 

k{c) = \ and U(c)=P + - -51u IN+scc . (10) 

1 — cp 1 — cp 

Then the derivative of 7r IN+scc (c) with respect to c is given by 

?4 + scc(c) = u IN+scc {k'(c)I + k(c)[I - cUicT^cUic))'} [I - cCZ(c)]- 1 , (11) 
where using ([TO)) after simple calculations we get 



Let us consider the point c = 0. Using pT|) . we obtain 

^n+scc(°) = -a(l -/3)u IN+scc +au IN+scc P. (12) 

One can see from the above equation that the PageRank of pages in IN+SCC with many 
incoming links will increase as c increases from zero, which explains the graphs presented in 

Next, let us analyze the total mass of the IN+SCC component. From (fT2|) we obtain 



IKn+scc(°)II = ~ a ( l - /3)u IN+scc + au IN+scc Fl = a(-l + (3+pi), 

where pi = u IN+scc Pl is the probability that a random walk on the hyperlink matrix stays 
in IN+SCC for one step if the initial distribution is uniform over IN+SCC. If 1 — (3 < pi 
then the derivative at is positive. Since dangling nodes typically constitute more than 
25% of the graph [19], and p\ is usually close to one, the condition 1 — < pi seems to be 
comfortably satisfied in Web samples. Thus, the total PageRank of the IN+SCC increases 
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in c when c is small. Note by the way that if (3 — then ||7r IN+scc (c)|| is strictly decreasing 
in c. Hence, surprisingly, the presence of dangling nodes qualitatively changes the behavior 
of the IN+SCC PageRank mass. 

Now let us consider the point c = 1. Again using (fTTj) , we obtain 

^IN + SCcll) = — ^ _ ^ U ™ + SCc[-^ — P ~ _ ^>S'l U IN + SCc] • (13) 

Note that the matrix in the square braces is close to singular. Denote by P the hyper-link 
matrix of IN+SCC when the outer links are neglected. Then, P is an irreducible stochastic 
matrix. Denote its stationary distribution by 7fi N +scc- Then we can apply Lemma I A. II from 
the singular perturbation theory to (fT3|) by taking 

A = P, eC = P — P — - ^S'lu IN+ sccj 

and noting that eCl = Rl + (1 — a — (3)(1 — /3) -1 ST. Combining all terms together and 
using 7f IN+scc l = ||7fi N+ scc|| = 1 and u IN+scc l = ||u IN+scc || = 1, from (|A.1[) we obtain 

II ' n N " a ' 

I I^IN + SCCV^ 



IN+SCCV-VII 1/3- P1 i ci ' 

H 71 iN + SCC- n '- L T IN + SCC*^ 

It is expected that the value of 7f IN+scc -Rl + 1 ^^ a 7f IN+scc 5 l l is typically small (indeed, in 
our dataset INRIA, the value is 0.022), and hence the mass 
l 7r iN+scc(c)|| decreases very fast as c approaches one. 

Having described the behavior of the PageRank mass ||tt in+scc (c)|| at the boundary 
points c = and c = 1, now we would like to show that there is at most one extremum 
on (0,1). It is sufficient to prove that if ||7r( N+scc (co)|| < for some Co € (0,1) then 
ll 7r i / N+scc( c )ll — for all c > Co. To this end, we apply the Sherman-Morrison formula to 
([9]), which yields 

_ ^u m+scc [I - cP]^Sl _ 

7TIN + SCC(C) = 7T IN+SCC (C) + c2a 7T IN+SCC (C), (14) 

+ l — c/3 u in + scc^ ~ C -T J <5J- 

where 

(1 — c)a 

ttin+scc(c) = — ^-u IN+scc [J - cP}- 1 . (15) 

1 — cp 

represents the main term in the right-hand side of (|14|L (The second summand in (|14|) is 
about 10% of the total sum for the INRIA dataset for c = 0.85.) Now the behavior of 
7r IN+scc (c) in Figure 3 can be explained by means of the next proposition (see Appendix IA. 21 
for the proof). 

Proposition 2 The term \ \ttw+scg(c) \ | given by \15\) has exactly one local maximum at some 
c G [0, 1]. Moreover, |K +scc (c)|| < for c € (c , 1]. 

We conclude that | |7r IN+scc (c)| | is decreasing and concave for c G [co, 1], where ||vr( N+scc (co)| | = 
0. This is exactly the behavior we observe in the experiments. The analysis and experiments 
suggest that c is definitely larger than 0.85 and actually is quite close to one. Thus, one 
may want to choose large c in order to maximize the PageRank mass of IN+SCC. However, 
in the next section we will indicate important drawbacks of this choice. 
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5 PageRank mass of ESCC 



Let us now consider the PageRank mass of the Extended SCC component (ESCC) described 
in Section [U as a function of c E [0,1]. Subdividing the PageRank vector in the blocks 
7r = [7r Purc0 uT tTbscc], from ([5]) we obtain 

lk ES cc(c)|| - (1 - c) 7 u ESCC [7 - cT^l, (16) 

where T represents the transition probabilitites inside the ESCC block, 7 = \ESCC\/n, and 
u ESC c is a uniform probability row- vector over ESCC. Clearly, we have ||7r ESCO (0)|| = 7 and 
||ti"escc (1)1 1 = 0. Furthermore, by taking derivatives we easily show that ||7r ESCC (c)|| is a 
concave decreasing function. In the next proposition (proved in the Appendix) , we derive a 
series of bounds for ||7r ESCC (c)||. 

Proposition 3 Let \\ be the Perron- Frobenius eigenvalue of T , and let p\ = u ESCC Tl be 
the probability that the random walk started from a randomly chosen state in ESCC, stays 
in ESCC for one step. 

(i) If pi < Ai then 

\\*escc(c)\\ < ^— A CG(0,1). (17) 
1 — cAi 

(ii) If 1/(1 - Pl ) < u ESCC [I - T}- 1 1 then 

\\^sac(c)\\ > ?p—°L, ce(0,l). (18) 



The condition p\ < Ai has a clear intuitive interpretation. Let 7r E scc be the probability- 
normed left Perron- Frobenius eigenvector of T . Then 7r ESCC , also known as a quasi- stationary 
distribution of T, is the limiting probability distribution of the Markov chain given that the 
random walk never leaves the block T (see e.g. 20 ). Since 7r ESCC T = Ai, the condition 
Pi < Ai means that the chance to stay in ESCC for one step in the quasi-stationary regime 
is higher than starting from the uniform distribution u ESCC . Although p\ < Ai does not 
hold in general, one may expect that it should hold for transition matrices describing large 
entangled graphs since quasi-stationary distribution should favor states, from which the 
chance to leave ESCC is lower. 

Both conditions of Proposition [3] are satisfied in our experiments. With the help of the 
derived bounds we conclude that ||7r ESCC (c)|| decreases very slowly for small and moderate 
values of c, and it decreases extremely fast when c becomes close to 1. This typical behavior 
is clearly seen in Figure [U where ||7T ESCC (c)|| is plotted with a solid line. The bounds 
are plotted in Figure Q] with dashed lines. For the INRIA dataset we have p\ = 0.97557, 
Ai = 0.99954, and for the FrMathlnfo dataset we have p 1 = 0.99659, A x = 0.99937. 

From the above we conclude that the PageRank mass of ESCC is smaller than 7 for any 
value c > 0. On contrary, the PageRank mass of Pure OUT increases in c beyond its "fair 
share" 8 = \PureOUT\/n. With c = 0.85, the PageRank mass of the Pure OUT component 
in the INRIA dataset is equal to 1.95(5. In the FrMathlnfo dataset, the unfairness is even 
more pronounced: the PageRank mass of the Pure OUT component is equal to 3. 445. This 
gives users an incentive to create dead-ends: groups of pages that link only to each other. 
Clearly, this can be mitigated by choosing a smaller damping factor. Below we propose one 
way to determine an "optimal" value of c. 
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Figure 4: PageRank mass of ESCC and bounds, INRIA (left) and FrMathlnfo (right) 



Let v be some probability vector over ESCC. We would like to choose c = c* that satisfies 
the condition 

||7r ESC c(c)|| = ||vT||, (19) 

that is, starting from v, the probability mass preserved in ESCC after one step should 
be equal to the PageRank of ESCC. One can think for instance of the following three 
reasonable choices of v: 1) ttt, the quasi-stationary distribution of T, 2) the uniform vector 
Uescc , and 3) the normalized PageRank vector 7r ESCC (c) / 1 |7r ESCC (c) 1 1 . The first choice reflects 
the proximity of T to a stochastic matrix. The second choice is inspired by definition of 
PageRank (restart from uniform distribution), and the third choice combines both these 
features. 

If conditions of Proposition [3] are satisfied, then (fl~7|) and (fl~8|) hold, and thus the value 
of c* satisfying (|19p must be in the interval (ci, C2), where 

(1 - ci)/(l -pici) = ||vT||, (1 - c 2 )/(l - Aic 2 ) - ||vT||. 

Numerical results for all three choices of v are presented in Table [1] 



V 


c 


INRIA 


FrMathlnfo 


T^ESCC 


Cl 




0.0184 


0.1956 




C2 




0.5001 


0.5002 




C* 




.02 


.16 


U B scc 


Cl 




0.5062 


0.5009 




C2 




0.9820 


0.8051 




C* 




.604 


.535 


tTescc/I |tTescc 1 


V(l + 


Ai) 


0.5001 


0.5002 




1/(1 + 


Pi) 


0.5062 


0.5009 



Table 1: Values of c* with bounds. 

If 

v — ^escc then we have ||vT|| — Ai, which implies C\ = (1 — Ai)/(1 — A1P1) and 
c 2 = l/(Ai + 1). In this case, the upper bound c 2 is only slightly larger than 1/2 and c* 
is close to zero in our data sets (see Tabel [TJ . Such small c however leads to ranking that 
takes into account only local information about the Web graph (see e.g. [21]). The choice 
v = 7T ES cc does not seem to represent the dynamics of the system; probably because the 
"easily bored surfer" random walk that is used in PageRank computations never follows a 
quasi-stationary distribution since it often restarts itself from the uniform probability vector. 



For the uniform vector v = u ESCC , we have ||vT|| = p±, which gives Ci,C2, c* presented in 
Table [TJ We have obtained a higher upper bound but the values of c* are still much smaller 
than 0.85. 

Finally, for the normalized PageRank vector v = 7r ESCC /||7r ESCC ||, using (| 16[> . we rewrite 
(fl9]l as 

7 7 2 (1 — c ) —i 

||TESCo(c)|| = J, ——n ESCC (c)Tl = - r^77U IN+scc [/ - CT] Tl, 

IFescc(c)|| ||7Tescc(c)|| 

Multiplying by ||7r ESCC (c)||, after some algebra we obtain 

||Wo(c)|| 2 =i|K scc (c)||-i^l. 

Solving the quadratic equation for || (c) 1 1 , we get 



II Mil M J T if c<l/2, 
||7r M oo(c)|| = r(c) = | ^ ifc>1/2 



Hence, the value c* solving (fT9|) corresponds to the point where the graphs of ||7r ES oo(c)|| and 
r(c) cross each other. There is only one such point on (0,1), and since ||7r ESCC (c)|| decreases 
very slowly unless c is close to one, whereas r(c) decreases relatively fast for c > 1/2, we 
expect that c* is only slightly larger than 1/2. Under conditions of Proposition r(c) first 
crosses the line 7(1 — c)/(l — Aic), then ||7Pr(c)||i, and then 7(1 — c)/(l — Pic). Thus, we 
yield (1 + A^ -1 < c* < (1 + p 1 )~ 1 . Since both Ai and pi are large, this suggests that c 
should be chosen around 1/2. This is also reflected in TabelQ] 

Last but not least, to support our theoretical argument about the undeserved high rank- 
ing of pages from Pure OUT, we carry out the following experiment. In the INRIA dataset 
we have chosen an absorbing component in Pure OUT consisting just of two nodes. We have 
added an artificial link from one of these nodes to a node in the Giant SCC and recomputed 
the PageRank. In Table [5] in the column "PR rank w/o link" we give a ranking of a page 
according to the PageRank value computed before the addition of the artificial link and in 
the column "PR rank with link" we give a ranking of a page according to the PageRank 
value computed after the addition of the artificial link. We have also analyzed the log file of 
the site INRIA Sophia Antipolis (www-sop . inria.fr) and ranked the pages according to the 
number of clicks for the period of one year up to May 2007. We note that since we have the 
access only to the log file of the INRIA Sophia Antipolis site, we use the PageRank ranking 
also only for the pages from the INRIA Sophia Antipolis site. For instance, for c = 0.85, the 
ranking of Page A without an artificial link is 731 (this means that 731 pages are ranked 
better than Page A among the pages of INRIA Sophia Antipolis). However, its ranking 
according to the number of clicks is much lower, 2588. This confirms our conjecture that 
the nodes in Pure OUT obtain unjustifiably high ranking. Next we note that the addition of 
an artificial link significantly diminishes the ranking. In fact, it brings it close to the ranking 
provided by the number of clicks. Finally, we draw the attention of the reader to the fact 
that choosing c = 1/2 also significantly reduces the gap between the ranking by PageRank 
and the ranking by the number of clicks. 

To summarize, our results indicate that with c = 0.85, the Pure OUT component receives 
an unfairly large share of the PageRank mass. Remarkably, in order to satisfy any of the three 
intuitive criteria of fairness presented above, the value of c should be drastically reduced. 
The experiment with the log files confirms the same. Of course, a drastic reduction of c also 
considerably accelerates the computation of PageRank by numerical methods [22, 5, 23 . 
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c 


PR rank w/o link 


PR rank with link 


rank by no. of clicks 


Node A 








0.5 


1648 


2307 


2588 


0.85 


731 


2101 


2588 


0.95 


226 


2116 


2588 


Node B 








0.5 


1648 


4009 


3649 


0.85 


731 


3279 


3649 


0.95 


226 


3563 


3649 



Tabic 2: Comparison between PR and click based rankings. 
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Appendix 

A.l Results from Singular Perturbation Theory 

Lemma A.l Let A(e) — A — eC be a perturbation of irreducible stochastic matrix A such 
that A(e) is sub stochastic. Then, for sufficiently small e the following Laurent series expan- 
sion holds 

[I - A(e)]- 1 = -X_! + X + eX 1 + ... , 

£ 



with 
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where [i is the stationary distribution of A. It follows that 



[I - A(e)]~ 1 = 



fjLeCl 



1/i + 0(1) as e -► 0. 



(A.l) 



Lemma A. 2 Let A(e) = A + sC be a transition matrix of perturbed Markov chain. 

The perturbed Markov chain is assumed to be ergodic for sufficiently small e different 
from zero. Let the unperturbed Markov chain (e = 0) have m ergodic classes. Namely, the 
transition matrix A can be written in the form 



A = 



Ax 


Li 







A,; 

L u 




E 



Then, the stationary distribution of the perturbed Markov chain has a limit 

lim n(e) = [viHi ■ ■ ■ v m \i m 0], 

where zeros correspond to the set of transient states in the unperturbed Markov chain, fii is 
a stationary distribution of the unperturbed Markov chain corresponding to the i-th ergodic 
set, and Ui is the i-th element of the aggregated stationary distribution vector that can be 
found by solution 

vD = v, v\ = 1, 
where D = MCB is the generator of the aggregated Markov chain and 



M 











fJ-m 



B 







with 4>i = [I - E]~ 1 L. l l- 

The proof of this lemma can be found in [TTJ [TJ1 fTJ] . 



A. 2 Proofs 

Derivation of {Pp. First, we observe that if 7r IN+scc and 7r DN l are known then it is straight- 
forward to calculate 7r OU T- Namely, we have 



1 



tt dn 1- )l T [I-cQ}- 1 . 



: cR[I-cQY 



Therefore, let us solve the equations (J7J) and ©. Towards this goal, we sum the elements 
of the vector equation |(5J), which corresponds to the postmultiplication of equation ([8]) by 
vector 1. 



1 



TlN+SCcCST + 7T DN 1 — 7T DN 11 1 — 

n n 

Now, denote by nour, nscc and n^N the number of pages in OUT component, SCC 
component and the number of dangling nodes. Since 1 T 1 — n^^, we have 



n — cudn 



,cSl + 



1 



: n DN ). 
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Substituting the above expression for 7r DN l into ((7j), we get 



i-cP sir 

n — cnrjN 



1 — C _ rp 1 



n — criDN n 



-udnI 



which directly implies ©■ 

Proof of Proposition [7] First, we note that if we make a change of variables e = 1 — c 
the Google matrix becomes a transition matrix of a singularly perturbed Markov chain as 
in Lemma [A . 2 1 with C = -^11 T — P. Let us calculate the aggregated generator matrix D: 

D = MCQ = -11 T Q - MPQ. 

n 

Using MP — M, MQ = I, and Ml = 1 where vectors 1 are of appropriate dimensions, we 
obtain 

D=-U T Q-I=-l[n 1 + l[I-T]- 1 R 1 l,--- , n m + 1[I - f ]" 1 i? m l] - I, 
n n 

where rij be the number of nodes in the block Qi, i — 1, ...,m. Since the aggregated 
transition matrix D + I has identical rows, its stationary distribution v is just equal to these 
rows. Thus, invoking Lemma I A. 2 1 we obtain ([4]). 

Proof of Proposition [H Multiplying both sides of (fl~5|) by 1 and taking the derivatives, 
after some tedious algebra we obtain 

Pi'n + scc(c)II = - o(c) + llw(c)||, (A.2) 

where the real- valued function a(c) is given by 

a(c) = -^_u IN+scc [/ - cP]- 1 ^ -P][I- cP]~ l l. 
1 — cp 

Differentiating (|A.2[) and substituting 1 _ e J |7r,sc'C'(c)|| from (I A.2|) in the resulting expression, 
we get 

IK+sccMII = |-a'(c) + + ||^ CC (c)||. 

Note that the term in the curly braces is negative by definition of a (c). Hence, if ||7i' I ' N+scc (c)|| < 
for some c E [0, 1] then | I^n+sccC ) 1 1 < f° r this value of c. 

Proof of Proposition fj| (i) The function /(c) = 7(1 — c)/(l — Aic) is decreasing and 
concave, and so is ||7r ESCC (c)||. Also, ||7r ESCC (0)|| = /(0) = 7, and ||7r BSOO (l)|| = /(l) = 0. 
Thus, for c G (0, 1), the plot of ||7T ESCC (c)|| is either entirely above or entirely below /(c). In 
particular, if the first derivatives satisfy | |7r ESCC (0)| | < /'(0), then ||7r ESCC (c)|| < /(c) for any 
c G (0, 1). Since /'(0) = 7(Ai — 1) and ||7r ESCC (0)|| = 7(pi — 1), we see that p\ < Ai implies 

The proof of (ii) is similar. We consider a concave decreasing function g(c) = 7(1 — 
c)/(l — pic) and note that g(Q) = 7, g(l) = 0. Now, if the condition in (ii) holds then 

ff'(l) > IkMoaWII. which im P lies d- 
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